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1 La justification du choix du sujet : 


Conscient de la prépondérance des énergies renouvelables qui constituent une 
solution miraculeuse de la problématique mondiale a savoir le réchauffement climatique et en 
programmant de faire la filière « environnement » dans la troisième année de ma formation 
ingénieur au sein de l’EMI .Pour ces deux simples raisons, j’opte pour ce sujet qui étudie 
l'optimisation du système d’une turbine éolienne tout en désirant d'enrichir mon bagage 
consternant tout ce qui est écologique. 


2 Application au Maroc: 


Il y'a quelques années, en 2014 plus précisément, le Maroc importait environ 90 Oé 
de ses besoins énergétiques. Le pays avait alors pris l'engagement d'améliorer la sécurité de 
l’approvisionnement en électricité, en réduisant la dépendance aux exportations d'énergie 
électrique et en augmentant l’utilisation des sources renouvelables. Depuis, la production 
d'énergie solaire et éolienne terrestre du Maroc s’est nettement développée. Fin 2019, la 
capacité d'énergies renouvelables du Maroc a atteint 3 685 MW, dont 1 215 MW d'énergie 
éolienne. Onze centrales éoliennes contribuent à cette nouvelle production d'énergies 
renouvelables, la centrale de Tarfaya dispose de la plus grande capacité, avec 301 MW. 


3 Des perspectives pour un élargissement de ce projet: 


Une éolienne en mer ou éolienne offshore est une éolienne implantée au large des côtes 
plutôt que dans les terres, pour mieux utiliser l'énergie du vent et produire de l'électricité grâce à une 
turbine et à un générateur électrique. Il existe deux principaux types d'éoliennes en mer : les éoliennes 
fixes, qui sont implantées sur des hauts-fonds, et les éoliennes flottantes qui offrent l'avantage de 
pouvoir être construites sur terre et implantées dans des zones où la profondeur des fonds marins ne 
permet pas la construction de fondations. Les coûts de production et d'installation sont supérieurs à 
ceux des éoliennes terrestres classiques, mais leur facteur de charge est plus élevé car elles bénéficient 
d'un vent plus fort et plus régulier. En 2019, Plus de 5 000 éoliennes sont connectées, en Europe, dont 
502 ont été installées dans l'année. La puissance crête cumulée représente 22,1 GW, répartis sur 110 
parcs dans douze pays européens, Royaume-Uni (44 %) et Allemagne (34 %) en tête. Ces deux pays, ainsi 
que le Danemark, les Pays-Bas et la Belgique, abritent en 2017 98 % du parc européen, favorisé par 
une me du Nord peu profonde et régulièrement bien ventée. 


4 Systèmes de Conversion Eoliens : 
4.1 Introduction : 


Le vent est une source d'énergie renouvelable, économique, exploitable avec un bon niveau de 
sécurité et respectueuse de l’environnement. Dans le monde entier, les ressources d'énergie 
éolienne sont pratiquement illimitées. Les récents développements technologiques dans les 
domaines des turbines éoliennes à vitesse variable, en électronique de puissance et en commande 
de machines électriques tendent à rendre l'énergie éolienne aussi compétitive que l'énergie 
d'origine fossile. 

4.2 Classement des Turbines Éoliennes : 


La technologie qui permet de profiter de l'énergie du vent a évolué sous diverses 
formes et types de machines. La structure de base des turbines éoliennes consiste aujourd’hui 
en un rotor pour capter l'énergie du vent en la transformant en énergie en rotation, un 
système d’engrenage pour démultiplier la vitesse de rotation du rotor, une machine électrique 
pour convertir l'énergie mécanique en électricité. Il existe différentes façons de classer les 
turbines éoliennes mais celles-ci appartiennent principalement à deux groupes selon 
l'orientation de leur axe de rotation :celles à axe horizontal et celles à axe vertical. 
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4.2.1 Turbines Eoliennes à Axe Horizontal (HAWT) : 


Une turbine à axe de rotation horizontal demeure face au vent, comme les hélices des avions 
et des moulins à vent. Elle est fixée au sommet d'une tour, ce qui 
lui permet de capter une quantité plus importante d'énergie éolienne. La plupart des 


éoliennes installées sont à axe horizontal. Selon son nombre de pales, une HAWT est dite 
mono-pale, bipale, tripale ou multi-pale.Une éolienne mono-pale est moins coûteuse car les 
matériaux sont en moindre quantité 
et, par ailleurs, les pertes 
aérodynamiques par poussée (drag) 
sont minimales. Suivant leur 
orientation en fonction du vent, les 
HAWT sont dites en « amont » (up- 
wind) ou en « aval » (down-wind). Les 
premières ont le rotor face au vent ; 
puisque le flux d'air atteint le rotor 
sans obstacle, le problème de « 
l'ombre de la tour » (tower shadow) 
est bien moindre. Néanmoins, un 
mécanisme d'orientation est 





essentiel pour maintenir en 
permanence le rotor face au vent. 


4.2.2 Turbines Eoliennes à Axe Vertical (VAWT) : 


L’axe de rotation d’une VAWT est vertical par rapport au sol et perpendiculaire à la 
direction du vent. Ce type de turbine peut recevoir le vent de n'importe quelle direction,ce qui 
rend inutile tout dispositif d'orientation. Le générateur et la boite d’engrenages sont disposés 
au niveau du sol, ce qui est plus simple et donc économique. La maintenance du système est 
également simplifiée dans la mesure où elle sefait au sol. Ces turbines ne disposent pas de 
commande d'angle de pale comme certaines HAWT. 
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5 Optimisation d’un Système de Conversion Eolien : 
5.1 Modèle du Système : 








Le système étudié est composé d’une turbine éolienne à axe horizontal 
tripale qui prend l'énergie de la masse d’air en mouvement, d’une boite de vitesses élévatrice 
qui adapte les vitesses de rotation de l’éolienne et du générateur, d’une machine synchrone à 
aimants permanents pour la conversion électromécanique, d’un pont à diodes qui fait la 
conversion électrique AC/DC et d’un groupe de batteries pour le stockage d’énergie. La charge 
est supposée consommer toute l'énergie produite. 








ns Système Mécanique: 


La puissance mécanique Pt qu’une turbine éolienne peut extraire d’une 
masse d'air traversant la surface balayée par son rotor est: 


E RE Ay 


p est la densité de fair (Kg/m°), A est la surface balayée par de rotor de l’éolienne (m2), v est la vitesse du vent et 
C, est le coefficient de puissance de la turbine. Ce dernier dépend du rapport de vitesses À (ou TSR, tip speed 
ratio). 


QR 
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TSR =A= 





La caractéristique non linéaire du coefficient de puissance C, peut s’approximer par une fonction rationnelle. La 
forme rationnelle proposée dans l'équation a l'avantage de montrer de façon explicite des informations telles 
que le TSR maximal pour un C, positif, ào et la valeur approximative du TSR optimal pour Cp maximal A = (Ao-a). 
Une simple régression de moindres carrés peut s'utiliser pour ajuster les coefficients G et a. 


GA. À) 


OT a + (4-2) 


Pour adapter la vitesse de rotation relativement lente de la turbine éolienne à celle du générateur, une boite 
d’engrenage (boite de vitesses) peut s'utiliser. Ce système de transmission mécanique dans laquelle M représente 
le rapport de transformation (ou transmission) de la boite, Q est la vitesse de rotation de l'arbre lent de la 
turbine éolienne et ( celle dela machine électrique. 


Q.=M.Q 


Apres le calcul on obtient : 


Q (A v-R 0) 


(a v) +A v- R Q) b 





5.3 Système Electrique: 


Le générateur est une machine synchrone à aimants permanents qui est modélisée simplement par 
une source de tension avec une impédance en série. Le circuit équivalent et le diagramme de Behn- 
Eschenburg sont montrés au dessous. Les composantes fondamentales pour la tension u, et le courant i, 
sont supposées en phase car la charge est un simple redresseur à diodes. 
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Les relations découlant de ce modèle simplifié de la machine sont les suivantes : 


E: valeur efficace de la composante 
fondamentale de tension induite parles aimants dans le bobinage du stator de la machine (f.e.m.). 
r uo di W ;: flux crête reçu par une bobine du stator 
= ir = roi venant des aimants. 
V2 J2 2 OO vitesse de rotation du champ 
magnétique (pulsation électrique ;œ = 27 f). 








ech o > E= 


w= p üg. (Qc = M Q, p: nombre de paires de pôles de la machine. 


QG: vitesse de rotation de l’arbre du générateur 
(œ= p Q6). 


DE E EE 


h; 


— 


Q: vitesse de rotation de l'arbre de la turbine (Qs =M Q). 
M: rapport de la boite de vitesses (multiplicatrice ou élévatrice). 


Grâce à quelques opérations algébriques sur le système précédent, on obtient une expression de la puissance de 
la machine définie par les paramètres R; et Let par la tension de batterie u, qui est une grandeur fixe dans ce cas. 





5.4 Paramètres du Système: 


Les caractéristiques mécaniques de la turbine éolienne, les paramètres de la fonction 
d’approximation du coefficient de puissance et les valeurs nominales et les paramètres caractéristiques du 


générateur à aimants permanents sont résumés dans Le tableaux 2.1, 2.2 et 2.3 respectivement : 


Tableau 2.1. Paramètres de la turbine éolienne Tableau 2 2 Coeficients de la fonction d’ approximation dl CI 








Paramètre Valeur 

Rayon (À) 1.8 m Deen in al 
Surface de balayage (4) 10.18 m° dume VW 
Coefficient de puissance maximal (LC ais) 0.42 Gan (C) UI 
TSR optimal LU" 68 Wee e 
Vitesse du vent nominale (vy) 12 m's Facteur D Lp 
Vitesse de rotation nominale (Qy 700 tr/mn TSR max () di 





Tableau 2.3. Générateur à aimants permanents 





Paramètre Valeur 8 
Couple nominal (Ty) S Nm 

Vitesse de rotation nomunale (Lly) 2000 tr/mn (210 rad's) 
Puissance nonunale(Py) 1680 W (2.23 HP) 
Tension nominale (vy) 110 Vag 

Resistance du bobinage de stator (Re) 0.9585 Q 

Inductance de bobinage de stator (Lẹ) 5.25 mH 

Flux induit par les aimants (#) 0.1827 Wb 

Nombre de pairs de poles (p) — 





La figure suivante montre comment la puissance évolue en fonction de la vitesse de 
rotation de la machine, avec plusieurs valeurs pour la tension de la batterie et une 
valeurde M constante. On peut observer sur la figure ou avec des tensions de batterie 
faibles, la machine peut commencer à fournir de la puissance à une vitesse de rotation 
basse. Cependant, avec une tension de batterie réduite, la valeur maximale de 
puissance produite par le système est aussi plus faible. 





La figure au dessous montre l’évolution de la puissance de la machine en fonction de la vitesse de 
rotation pour plusieurs valeurs du rapport de transformation de la boite de vitesses M,avec une tension 
de batterie fixe. 
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On peut voir à partir de la figure que l'effet plus important relié à M est ou avec un rapport plus élevé, la 
machine commence à fournir de la puissance à une vitesse de rotation inférieure. Ceci est intéressant 
pour profiter d’une vitesse de vent de démarrage plus faible, pour élargir la plage de vitesses de vents du 
système. Cependant, en même temps, la vitesse à laquelle le système décroche à cause de la surcharge 
(P > Pom) est aussi plus faible, ce qui diminue la plage de vitesses du côté des valeurs supérieures. En 
conséquence, il est important de bien choisir la valeur de M de façon à maximiser la plage de vitesses du 
système : il doit être assez élevé pour faire démarrer le système à des vitesses faibles, mais assez réduit 
pour permettre au générateur d'atteindre les vitesses supérieures. 





= Dans la suite, un problème d’optimisation du système de génération éolien présenté 
est défini pour maximiser la puissance produite du système en cherchant les valeurs 
optimales du rapport de transformation de la boite de vitesses et de la tension de 
batterie. 


Kë Optimisation du système : 


Les équations de puissance mécanique et la puissance de générateur en régime permanent du système sont ici 
les expressions analytiques qui permettent la formulation de l'objectif principal du problème d’optimisation. Le 
point de fonctionnement permanent du système se trouve àl'intersection des deux courbes représentant ces 
deux puissances dépendent de la vitesse du vent, et de la valeur du rapport de transformation de la boite de 
vitesses M et de la tension de batterie u, 


Le problème d'optimisation peut alors être posé de la manière suivante. Il consiste à trouver le jeu de paramètres 
permettant au système éolien de maximiser la puissance produite sur la plage de vitesse du vent : 


max P 
[M vi! 


5.5.1 Contraintes: 


Les boites d’engrenages parallèles à un étage ont des rapports de transformation maximaux de 1:5.En 
associant plusieurs étages, le rapport de transmission soit plus élevé. Ce rapport doit respecter une contrainte 
imposée par les valeurs nominales des vitesses de rotation de la turbine et du générateur Qsenn et Qu. Un 
système de faible taille utilise une éolienne qui tourne relativement vite, la vitesse maximale d’un générateur 
électrique de faible puissance est de 3600 t/m. Le rapport entre la vitesse du générateur et celle de la turbine 
OGenn /Q, risque donc d’être plus faible que le rapport maximal envisageable. Cette valeur devient une borne 
supérieure pour M. 


Les valeurs nominales du générateur imposent les limites de tension et de courant. Il est supposé que ces 
restrictions sont suffisantes pour maintenir la puissance générée en- dessous la puissance nominale et que la 
turbine éolienne peut délivrer toute la puissancemécanique pour les vitesses de vent faibles et modérées (v < vy). 
Quand le vent atteint la vitesse maximale Ivan), la petite éolienne s’auto protège des vents destructeurs en 
sortant dela direction du vent. 


L’équation qui modélise la puissance du générateur n’est valide qu’à partir du moment où la tension 
induite est supérieure à la tension seuil imposée par la tension de labatterie pour que les diodes du pont 
soient passantes. Cette condition impose une vitessede rotation minimale pour que le générateur commence à 
fournir de la puissance à la charge. La tension de batterie oblige indirectement à une vitesse de ventminimale 
(Vout-in): 
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En conséquence, la formalisation du problème d'optimisation proposé est la suivante : 
Trouver les paramètres M et û, tels que : 


max P 
[M u,] 


Avec les contraintes : 





5.5.2 Problème d'optimisation: 


La recherche analytique de la solution du problème ainsi défini pose néanmoins quelques 
difficultés : 
1) La réduction à une seule équation n’est pas possible. 


L’exploitation de l’équation d’égalité des puissances ne permet pas d’extraire la seule 
variable indépendante qu’elles ont en commun, la fréquence de fonctionnement (œ). De 
ce fait, 1l n’est pas possible d’obtenir une expression de la puissance à maximiser à 
partir des seuls paramètres d’optimisation. 


2) La paramétrisation avec la seule variable indépendante non contrôlable (v) ne 
mène pas à une solution unique. 


Pour une valeur de la vitesse de vent donnée, 1l y a une vitesse de la turbine qui 
correspond à une production maximale de puissance éolienne, cette vitesse est notée Q“. 
Pour chaque valeur du rapport de transformation de vitesse, M, il y correspondra une 
fréquence de fonctionnement du générateur électrique notée 

œ* donnée par : 





5.5.3 Adaptation du Problème d’Optimisation: 


Pour remedier a ce problme ,on utilise cette fois seulement la fonction de la puissance mécanique, une forme plus 
simple est utilisée. Pour réduire les degrés de liberté du système, la vitesse du vent, seule variable non 
contrôlable du système, est sortie de la formulation mathématique par l’utilisation d’une forme optimale. 





Si le rapport de vitesse À est maintenu à sa valeur optimale À, le coefficient de puissance est toujours à sa valeur 
maximale Cpm = C.L Donc, la puissance del’éolienne est aussi à sa valeur maximale. On obtient donc une forme 
analytique de la puissance mécanique maximale de laturbine éolienne en fonction de sa vitesse de rotation Q 

uniquement. 





La surface entre les courbes de puissance mécanique idéale et la puissance produite parla machine est : 


OD: 
4=4-4,= [{PQ)-P,(&.M,Q)}dQ 


L'objectif du nouveau problème d'optimisation est de rapprocher les deux courbes ; 
donc de minimiser la différence entre ses aires : 





La puissance mécanique idéale de la turbine éolienne P; varie selon la vitesse et atteintsa valeur nominale Py à 
la vitesse de vent nominale vy . Il y a cependant un rang de vitesses de vent entre v, et la valeur maximale (cut- 
off), où la puissance de la turbine éolienne doit être régulée de façon à ne pas dépasser Py. Pour les petites 
éoliennes, ceci est fait par le système de régulation aérodynamique de type stall. Qç est la vitesse de rotation de 
la turbine à laquelle la puissance arrive à Py. Il est important de noter que les vitesses Qc et Qp (vitesse nominale 
de rotation de l’éolienne)ne sont généralement pas égales (Qc < Qy). 
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l'intégration de P; donne une valeur fixe, qui dépend uniquement des caractéristiquesde l’éolienne. 
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Les restrictions physiques du système et les contraintes mathématiques de l'équation de la machine 
permettent de définir les limites d'intégration. La limite supérieure est obtenue soit par la valeur 
nominale de la vitesse de rotation de l’éolienne, soit par la valeur de la vitesse où la puissance électrique 
de la machine est supérieure à la puissance mécanique idéale ou à la puissance nominale. La limite 
inférieure est obtenue soit de la valeur minimale de fonctionnement du système, soit de la condition de 
positivité pour l'équation de la puissance, soit de la condition de puissance non imaginaire: 


Clinas = MIN Io < Liv: PaM, H, LI) SP: P mM, H, (à) < PN } 





Omin = Max {Q > Omin sys ; PM Q-J(R, YY +E -|(pM YQF -G|-R,u >0; 


LG vert -| (pM YQF -u |>0} 


Les limites pour la tension du système et du rapport de transformation sont les mêmesque pour le 
problème précédent : 
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Une dernière contrainte utilisée est de limiter la puissance de la machine à P; pour éviterun surdimensionnement 
de la machine. Ceci a été fait pour toute la plage de vitesses de fonctionnement du système : 


PA M. Hs, (2) = : PQ) vie ISB LA 


Pour résoudre le problème d’optimisation précédent, une méthode de Monte Carlo a étéutilisée. 


La procédure de solution est : 


l. Choisir un nombre SP de paires (us, M) initiales dans les limites de l’espace retenu 

2. Vérifier les conditions pour les paires choisies et garder uniquement les paires qui 
satisfont les contraintes du problème (solutions faisables). 

3. Créer une fenêtre de recherche avec les valeurs minimales et maximales des 
solutions faisables trouvées [Umin Mi, Una Mmaxl. 

4, Choisir un vecteur de recherche r = fr,, rm] aléatoire ; chaque composant a une 
valeur entre 0 et 1 et estimer les variables d'optimisation par 

Hs Oli wu. u 


N Har miti min 


M| |0 n| |Ma Mul |M 


faihl 


Evaluer la faisabilité de la paire choisie et en cas favorable, 
Calculer l'intégrale Am numériquement pour chaque paire faisable. 


Garder les valeurs de u, M et A. 


EE A 


Répéter les étapes 4 à 8 un nombre de fois N avec un nouveau r à chaque 
itération. 
9. Arranger les N résultats antérieurs en ordre croissant. 


10. Garder les E premiers (meilleurs) résultats pour refaire une nouvelle fenêtre de 
recherche et répéter G fois les points 3 à 10. 


Les paramètres SP, N, et E sont des valeurs arbitraires. Ainsi, à la fin de la dernière itération de la procédure, la 
solution du problème se trouve à la première place des derniers résultats rangés. 


5.5.4 Résultats: 
Les paramètres du système sont toujours les mêmes que ceux du cas précédent. Dans le tableau, les solutions 
obtenues pour 5 cas sont montrées. Un programme sur MATLAB fut préparé et utilisé pour rechercher les solutions. 
Le nombre initial de candidats fut SP = 20. Pour chaque essai de Monte-Carlo, N = 20 individus furent testés. Le 
nombre de fois que les essais furent répétés est G = 10. Les différentes solutions obtenues se doivent à l'existence de 
multiples points optimaux locaux et à la caractéristique aléatoire des points initiaux choisis par la méthode. 
Cependant, tous ces points sont dans une même zone, on peut donc dire que la tension optimale u se trouve entre 
30 V et 33 V et que le rapport optimal de transformation de la boite de vitesses M se trouve entre les valeurs 2.1 et 
2.5. 


Cas A / Ai He M 
I 0.1919 20 9033 2.5338 
Il 02075 31.5285 2.3975 
II 0.2303 32.3938 2 2885 
IV 0.2444 32.5845 2.2381 
V 0.2769 32.9100 2.1322 





Il est remarquable que le cas I donne le meilleur résultat : la surface relative A / A; est la plus petite des cas 
réalisés, qui peut être considéré comme le cas optimal ; donc, les valeurs optimales de la tension du 
système et du rapport de transformation de la boite devitesses sont est u, = 30 V et M = 2.5. 


6 Conclusion: 

Une méthode pour l'optimisation d’un système isolé de conversion éolien de petitetaille est 
présentée et étudiée. L’objectif est de maximiser la puissance produite par un système simple sans commande 
mécanique ni électronique. La méthode est basée sur unmodèle simple sans pertes dans la transmission mécanique 
avec lequel on obtient les expressions de la puissance mécanique de la turbine éolienne et électrique de la machine. 
Une procédure analytique permet de trouver, pour chaque vitesse de vent, les valeurs de M en fonction de la 
tension de batterie, afin de maximiser la puissance produite. Ainsi,il est nécessaire d’établir un critère qui permette 
de définir une paire unique M et u, pourle système. Un critère de minimisation de la surface entre les courbes d’une 
puissance idéale de référence et la puissance de la machine a été utilisé pour toute la plage de variation de vitesse 
du vent. Ainsi, une procédure de solution par la méthode de Monte Carlo a permis de trouver une zone de points 
optimaux qui permet de maximiser la puissance générée par le système de conversion éolien. 


7 Efficacite des projets d'optimisation effectués par un bureau d'étude: 


Généralement les projets d'optimisation servent a plusieurs choses, parmi ces choses on 
trouve ` améliorer l'efficacité, fournir un accès à des informations pertinentes et exactes, meilleure 
Adaptabilité, améliorer la qualité des résultats, voire même suivi de la Performance et 
Responsabilité. Par conséquent, si ces projets vont être réalisé dans un bureau d'etude, l'efficacité 
d'optimisation sera garantie ce qui va augmenter la rentabilité. 


8 References : 


httos://Wwww.fedenerg.ma/2020/05/29/energie-eolienne-la-banaue-mondiale-met-en- 





avant-le-grand-potentiel-du-maroc/ 


httos://fr.wikipedia.org/wiki/Liste des parcs %C3%A9oliens au Maroc 


httos://fr.wikibedia.org/wiki/%C3%89olienne en mer 


httos://Wwww.connaissancedesenergies.org/fiche-bedagogique/eoliennes-en-mer- 





offshore 
https://www.voutube.com/watch?v=J7K5KDiale4 
https://www.voutube.com/watch?v=Z0wtmej2zEM 
https://www.voutube.com/watch?v=YWwHvaTwQERM 
https://www.voutube.com/watch?v=x9Gtdaq9oVf4 
https://www.voutube.com/watch?v=DKyviAdz8qc 

__ https://tel.archives-ouvertes.fr/tel-00344978 
https://core.ac.uk/download/pdf/78383521.pdf 


__ https://dspace.univ- 


adrar.edu.dz/jspui/bitstream/123456789/3422/1/OptimisationX20d%27un%20systXC3%A 
8me.pdf 


httos://WwWww.instagram.com/p/CLevmTWKG6 /?igshid=luugOnno5nxga 


